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Abstract 

It was proposed that a double quantum dot can be used to be a detector of spin bias. Electron 
transport through a double quantum dot is investigated theoretically when a pure spin bias is 
applied on two conducting leads contacted to the quantum dot. It is found that the spin polarization 
in the left and right dots may be induced spontaneously while the intra-dot levels are located within 
the spin bias window and breaks the left-right symmetry of the two quantum dots. As a result, a 
large current emerges. For an open external circuit an charge bias instead of a charge current will be 
induced in equilibrium, which is believed to be measurable according to the current nanotechnology. 
This method may provide a practical and whole electrical approach to detect the spin bias (or the 
spin current) by measuring the charge bias or current in a double quantum dot. 

PACS numbers: 85.75.-d, 85.35.-p, 73.21.La, 73.23.-b 
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I. INTRODUCTION 



Discovery and application of giant magnetoresistance (GMR) in metallic thin films marks 
the beginning of a new era of spintronics.-i^ Since then, people begin to exploit electron spin 
to replace the role of electron charge in electronic devices. As a counterpart of charge cur- 
rent, spin current, in which spin-up and spin-down electrons move coherently in opposite 
directions, has been attracted extensive interests.^ Various methods were proposed to gen- 
erate spin current,- and to explore the characteristics of the spin transport. Over last few 
years, search of spin current has made a great of progresses. It has been generated and de- 
tected successfully by various means, such as the optical injection,-i^ the magnetic tunnelling 
injection,-!^ or the spin Hall effect.-'^^ All these experiments focus on the optical measure- 
ment of spin accumulation near the boundaries of sample or electric measurement of the 
scattering effect induced by the spin current via spin-orbital coupling. There are also some 
proposals to measure spin current or spin polarized current,— i^'^^'i^ e.g. to measure the 
spin torque while a spin current flowing through a ferromagnetic-nonmagnetic interface,— 
or to detect the induced electric field by the spin current.— In all these methods, it always 
involves the optical, magnetic materials or impurities, magnetic field, or spin-orbit interac- 
tion. Up to now, it is still a challenge to detect the spin current efficiently, which has become 
a bottleneck of the development of the spintronics. 

When a spin current flows through a device, there always exists a spin bias between 
the two terminals of the device. A spin bias means that the chemical potentials of the 
two terminals are spin-dependent (see Fig.l). The spin bias is regarded as the driving 
force behind the spin current. When the circuit is open, the spin current has to be zero. 
Consequently the spin bias usually induce spin accumulation in equilibrium. When the 
circuit is connected, a spin current circulates. The relation between spin bias and spin 
current is very similar with the relation between the charge bias and charge current. On 
the charge transport, people often detect the charge bias to replace the measurement of 
the charge current. Correspondingly we can also measure the spin bias instead of the spin 
current. In this paper, we propose an effective method to detect the spin bias. 

The present proposal is a whole electric measurement of spin bias by means of a double 
quantum dot (DQD). It does not involve any optical or magnetic means, and even the spin- 
orbit interaction. The spin bias can be detected by measuring the (charge) bias. The DQD 
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can be regarded as an artificial molecule, and the electron numbers in DQD can be controlled 
very well. In last two decades, the electron transport through the DQD device has been 
extensively investigated.— >i'' DQD has also been proposed as a qubit,— a device to detect 
various tunnelling rates and spin flip rate,— and so on. Here we propose that a DQD can 
be applied to measure the spin bias or spin current. 

Let us flrst describe the working mechanism of DQD as a detector of spin bias. Consider a 
DQD coupled into two conducting leads. Suppose a spin bias be applied between the left and 
right leads. Our task is to measure this spin bias experimentally. The spin bias is deflned as 
the spin-dependent chemical potentials of the two leads with Hl^ = —f^Li = —f^Rt = f^Ri = V 
(see Fig.l).-° Assume that the left-dot level ei is set at zero and the right-dot level is at 
—U, where U is the intra-dot electron-electron (e-e) Coulomb interaction. This particular 
level position is chosen to demonstrate the physics in our proposal, and is not necessary at 
all in a general case. The left dot has a spin- up electron because of fi^i; > > fJ'Li, while 
the right dot, because of fiRi > eu + U > jiR-^ > e/j, is occupied by a spin-down electron, and 
its spin-up level is consequently pushed away to the higher energy eR + U and is empty (see 
Fig.l). The spin- up electron can then tunnel from the left lead via the two dots to the right 
lead (see Fig. la). Oppositely the spin-down electron can hardly flow from the right lead to 
the left lead because of the Pauli exclusion principle and the occupancy of the spin-down 
level in the right dot (see Fig. lb). This breaks the symmetry of the motion of spin- up and 
spin-down electrons in a pure spin bias. As a result, a (charge) current circulates. This 
induced current can be measured experimentally, and consequently be applied to measure 
the spin bias. 

The paper is organized as follows. In Section II, the model for the DQD and the gen- 
eral formalism for nonequilibrium Keldysh Green's function method are presented. The 
spin-bias-induced charge current J and the electron occupation numbers in the DQD are 
calculated. In Section III, we take the numerical investigation. The spin- dependent charge 
stability diagram in terms of the spin bias is obtained. In Section IV, the induced charge bias 
in an open circuit is numerically studied. Finally, a brief summary is presented in Section 
V. 
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II. MODEL AND FORMULATION 



In this section, we present the model Hamiltonian of this DQD and the general formalism 
of Keldysh Green's function technique for electron transport through the DQD. The DQD 
device is modelled by the following Hamiltonian, 

a,k,a a,(T 

a a, a' 

+ ^ tattlf^^daa + ^ tcd^^dn^ + H.c. (1) 

a,k,cr cr 

where al^f.^ {daka) and (ij^^ {daa) are the creation (annihilation) operators of electron with 
spin cr(='f, I) in the lead a(= L, R) and the dot a , respectively. Each dot has a single energy 
level and an intra-dot e-e interaction f/j„. In addition, the inter-dot e-e interaction Uex is 
also included. We emphasize that the system does not break the spin SU(2) symmetry, and 
the hopping coefficients and tc are spin-independent. 

Following the transport theory of Keldysh Green's function,— the electron current J^o- 
with the spin a from the lead a flowing into the dot a and the occupation number of electron 
Uaa at the level a, a can be expressed as, 

Jaa = -Im J [2/..GL.(e) + Gl,(e)] (2) 

naa = {dijaa) = ~^ j '^'^^^ai^) (3) 

where = 27r |tap5(e — ^ak)- /acr(e) = l/{exp[(e — fiaa)/kBT] + 1} is the Fermi-Dirac 
distribution of electrons in the leads. Because of the spin bias in the two leads, the chemical 
potentials for spin- up and spin-down electrons are not equal. GJ^„^(e) and G'^^^{e) in Eqs. 
(2) and (3) are the standard retarded and the Keldysh Green's functions of the QDs, they 
are the Fourier transformation of G^aa (^) 5 where 

We first solve the Green's functions g^(e) of the isolated DQDs system (i.e. ta = tc = 0). 
Consider that the spin bias V is less than the intra-dot e-e interaction Uin and the two- 
electron co-tunneling events can be ignored. g^(e) are obtained from the equation of motion 
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technique:— 



A A-U, 



ex 



A — Uin A — Uin — Ue. 

and g\ji^ = fi'/fXcr = O5 where a = R for a = L and a = L ior a = R, a =1 for a =| and 
a =t for a =[, A = e — ea — [na]Uex + {na} = n^ — [na], and [ua] is the integer part of 
Ha- Ua = + Hal is the total occupation number of electron in the dot a. After solving 
gJJ.(e) of the isolated DQDs, Gaacr(^) ^aaai^) the whole system can be obtained from 
Dyson and Keldysh equations:— 

G;(e)= I ^^^'^ ^^^'^ I = g;(6) + g;(e)S;G;(e), (5) 
G<(e) = I ^^^'^ 1 = G;(e)S<(e)G:(e). (6) 

^RLa ^ RRo 

Here the bold face letters (G, g, and S) represent the 2x2 matrix, and the self-energies 
S:'<('e) are: 

S;(e) = ' I , (7) 

te -iVnl2 

^^rJrA^) 

Eqs. (3, 4, 5, and 6) can be solved self-consistently. The (charge) current through the DQD 
is given by 

J = e{J^ + Jli) = -e{jR^ + Jri). 

Finally it is worth pointing out that the present problem can be solved by other means, 
for example, the rate equation method.— 



III. SPIN-DEPENDENT CHARGE STABILITY DIAGRAM AND CHARGE 
CURRENT 

Before presenting numerical results, we emphasize that the spin bias we apply to the 
DQDs device is a pure symmetric one without a (charge) bias, i.e. /^li + /^lj. = A^iit + A^rj. = 
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0.— So if the spontaneously spin-polarized occupations are not induced in the DQD, the 
charge current J must be zero because of the symmetric behaviors for the motion of spin- 
up electron and the spin-down electron. For example, in the case of a single quantum dot 
instead of DQDs applied by the pure spin bias, there is no spin polarization in the dot and 
the current is always zero as the spin up-down symmetry is retained. So, in the following, 
we first investigate the stability diagram of spin polarization and the spin-dependent charge 
density in the DQD. 

Fig.2a and b present the spin polarizations An^ (An^ = — Uai) of the left and right 
dots versus the levels and e^?, and Fig. 2c presents the occupation number of electron 

+^/?/2.^^ It is found that these quantities are determined by the relative energy levels of 
ei and e^. The spin polarization An^, is indeed non-zero and even quite large (i.e. near ±1) 
in some specific regions. Let us analyze the spin-dependent charge stability diagram (see 
Fig. 3a), which gives spin-dependent occupation numbers of electron as a function of and 
eR. If without the spin bias {V = 0), there are four domains (0, 1), (1, 1), (0,2), and (1,2) 
in the stability diagram (see the thin dashed curves in Fig. 3a), with {n,m) representing 
n and m electrons in the left and right dot. This type of charge stability diagram has 
been observed experimentally,—'^ and is well established. While the spin bias V is turned 
on and the level or eR locates between —V and +V, in addition of the four old spin- 
unpolarized domains {n, m) with a shift V of their boundaries shift, there appears four new 
spin-polarized domains, which are denoted by (t,!), (0, i), (l,i), and (1,2). The notation 
(t, 1), for example, represents an electron of spin- up in the left dot and a spin-unpolarized 
electron in the right dot. 

This spin-dependent charge stability diagram in Fig. 3a can be obtained by calculating 
the electrochemical potentials of the DQD or by analyzing the level's position relative to the 
spin-dependent chemical potentials fiaa- Consider the isolated DQDs device with = Tr = 
tc = 0. (i) The domain (0,1): when the equivalent level (e^ = + Uex) of the left dot is 
higher than fi^ and fiLi, and the right-dot's level eR satisfies eR < e^, /i/j-f, /x/ji < e_R + Uin 
(see Fig. 3b), the right dot is occupied by a spin-unpolarized electron and the left dot is 
empty, (ii) The domain (t,l): while fiLi < < yULT ^nd eR < fiRi,fiRi < eR + Uin (see 
Fig. 3c), a spin-up electron occupies the left dot and a spin-unpolarized electron is in the 
right dot. (iii) The domain (0, j): if ii > fJ'Lh f^H and fiRi^ < eR + Uin < fJ'Ri (see Fig.3d), 
the left dot is empty. For the right dot, a spin-down electron occupies the level eR because 
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of eji,eji + Uin < Urii then the spin-up level of the right dot is pushed to eR + Um which 
is over /x/j^, and so it is empty. Similarly, the other five domains can also be obtained. In 
the case of the finite coupling case F^, Vr, tc 7^ 0, the spin-polarized domains slightly extend 
to the spin-unpolarized domains as illustrated in the thin dotted lines in Fig. 3a. Numerical 
results for the spin polarizations An^ (Fig. 2a and b) and the occupation numbers of electrons 
nL + nR/2 (Fig.2c) are in a good agreement with the charge stability diagram in Fig.3a. The 
eight domains, including four spin-unpolarized and four spin-polarized domains, are clearly 
visible. 

In an alternative way, the stability diagram of Fig. 3a can also be deduced from the total 
energy of the DQD system and the electrochemical potentials. When the isolated DQD is 
in the states of iV = {N^, N^^, Nj^-, N^i), where Naa = or 1 is the index of the electron 
occupation number in the intra-dot level aa, its total energy Ejr is 

Et{N) = NLeL + NRen + NLNnUex 

+ iNL^NLi + NR^Nni)U^n, (9) 

with Na = + Nai- Cousidcr the fact that the occupation number in the intra-dot level 
aa is mainly effected by the lead (i.e. electron reservoir) aa. The grand thermodynamic 
potential Q at the zero temperature is 

Q{N) = Et{N) - N^fiL^ - NLim 

-iV/jtm - NRifiRi- (10) 

In the present system, the electron occupation number can change with the levels and eR. 
This is a grand canonical ensemble. Then the stablest state is one whose grand thermody- 
namic potential Q has the minimal values, and can be found straightforwardly. For the sake 
of convenience and intuition, we introduce the electrochemical potentials fiQDaa, following 
Ref.— . [iQDaa of the level aa is well defined, for example, 

/iQDLT(iV) = Et{N) - Et{N^ - 1, Nli, Nr^, Nr^). (11) 

Then the stablest states are the maximal values of for which four ^QDaa{N) are less than 
the corresponding chemical potentials fXaa- If two states of N, e.g. N = (0,0,1,0) and 
(0, 0, 0, 1), satisfy the above four equations, they are assumed to have the same probability 
to exist. A detailed analysis of iiqoaa versus the parameters ei and eR leads to establish the 
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same charge stability diagram as shown in Fig.3a. In fact, the electrochemical potentials 
IJ'QDaa ^rc cqual to the equivalent levels in the preceding paragraph. For example, 

/^QDLt(1, 0,1,0) = //Qi5L|(l,0,0, 1) =/iQDLi(0, 1,1,0) 
= IJ'QDLiiO, 1, 0, 1) = ez, + Uex = CL- 

In particular, there are only four equivalent levels, which are less than the numbers of /iQDaa- 
So it is convenient and intuitive to use the equivalent levels to deduce the stability diagram. 

With the spin-polarizcd stability diagram in mind, we turn to calculate the (charge) 
current J induced by the spin bias. Fig. 2d shows the current J as a function of the levels 
cl and e^. The current becomes quite large when both the left and right dots are spin 
polarized in the case of —V < ii — €r + Uin < V. The physical origin of generation of 
the current has been explained in detail in the introduction and as shown in Fig.l. We 
can establish a relation between the charge current and the spin bias in the two leads. In 
this way, we can detect the spin bias V by measuring the current J. In the following we 
calculate the current for various parameters. Fig.4a shows the current J versus the spin bias 
V for the inter- dot interaction Uex — 5. While V = 0, J is zero exactly. With the increase 
of V from zero, the current J first increases, reaches at a maximum, and then drops. J 
keeps a relatively large value even if V is comparable with the e-e interaction energy Uin- 
The origin of the drop is that the spin-polarizations in two dots decay while the current 
flows through the DQDs at the large V. In the absence of the inter-dot e-e interaction Uex, 
i.e. Uex = 0, the current increases monotonously with the spin bias V (see Fig. 4b). In this 
case the current J and the spin bias V have a one-to-one correspondence. Therefore the 
spin bias V can be deduced straightforward from the measured current. Fig.4c shows the 
current J as a function of the right-dot's level cr. When cr + Uin departs over a few F^ 
(e.g. \eR + Uin — ^l\ > SFq), J becomes very small because the tunneling process in Fig. la 
suppresses quickly when cr + Uin is not in alignment with . On the other hand, the 
tunneling process in Fig. la occurs frequently and J becomes large when eR + f/j„ is located 
near ex,. However, when sr + Uin = ^l, J may drop slightly and a dip emerges in the curve 
of J-eR, because the spin-polarization Aria is suppressed at the point. Fig.4d displays the 
current J as a function of temperature T. Here J depends on the temperature T slightly, 
and is quite large when T <V. 
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IV. CHARGE BIAS IN AN OPEN CIRCUIT 



In the preceding section, we calculated the charge current through a DQD induced by 
a pure spin bias. In an open circuit, the situation will be changed. At the time that a 
spin bias is turned on, a charge current will circulate. For an open circuit, the extra charge 
will accumulate in the two leads until the system reaches at a balance. As a result an 
extra charge bias VJ= instead of a charge current will be generated while the charge current 
vanishes. In this case, combination of the the spin bias V and the induced charge bias Ve 
will give the spin-dependent chemical potentials fiaa in the two leads 

/iLT = +V + Ve, (12a) 

m = -v + Ve, (12b) 

m = -V- K, (12c) 

m = +V- K. (12d) 
The bias Ve can be determined by the condition of 

J = (13) 

in equilibrium for an open circuit. Figs. 5a and 5b gives the bias Ve and Ve/V versus the spin 
bias V in the presence and absence of the intra-dot Coulomb interaction Uex- \Ve\ and \ Ve/V\ 
increase monotonously with V regardless of the value of Uex ■ This is different from the curve 
of J-V, in which J drops down for a large V while Uex 7^ (see Fig. 4a). This illustrates 
that it is more efficient to measure the induced bias Ve than to measure the induced current 
J. Fig. 5c shows the bias Ve as a function of the level ept- The bias \Ve\ always has a large 
value (e.g. \Ve/V\ > 0.1), even if + f/j„ is far away from e^. Notice that the current J 
is relatively small when |e/j + Uin — clI > SFq, (see Fig. 4c). The transmission coefficient (or 
the conductance) is also very small in this region. Correspondingly, Ve in an open circuit 
is still large. Therefore the induced bias Ve can be measured in an more extensive region. 
Fig.5d gives the temperature T dependence of the bias Ve, which is almost independent of 
the temperature T. Finally, we emphasize that \Ve/V\ is usually larger than 0.1 regardless 
of the values of the parameters V, e^, e^, T, etc. In the current technology, the bias in 
the order of O.lnV is measurable in experiment.— Therefore, if the spin bias V, i.e. the 
difference of the spin-up and spin-down chemical potentials (/iL-f — ^Li)/e, reaches to InV, 
the induced bias in the present calculation is large enough to be measured in experiment. 
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V. CONCLUSIONS 



In summary, we investigated the electron transport driven by a spin bias or pure spin 
current through a non-magnetic DQD. Except for the spin-unpolarized domains, several 
spin-polarized domains are found in the stability diagram with respect to the energy levels 
of two quantum dots. When both of the left and right dots are spin polarization, a large 
charge current J can be induced by applying of the pure spin bias. In particular, in an open 
circuit, the charge bias is induced to balance the spin bias, and is measurable in an extensive 
range of the parameters. Physically, a pure spin bias may drive electrons with different spin 
in opposite direction. If the system possesses the left-right symmetry or parity and does not 
break the time reversal symmetry it will circulate a pure spin current (or spin accumulation 
in an open circuit). When the energy levels in the two dots are not equal, the left- right 
symmetry or parity of the system is broken. A spin bias and strong Coulomb interaction 
can produce two spin polarized states in the two dots as we discussed in the spin-polarized 
charge stability diagram. As a result the currents with different spins in opposite directions 
will not be equal any more. Consequently this pure spin bias generates a charge current 
through the DQD. This property may provide a practical approach to detect the spin bias 
in DQD by measuring the charge bias or charge current. 
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FIG. 1: (color online) (a) ((b)) The sehematic plots illustrate an spin- up (spin-down) electron 
tunneling from the left (right) to the right (left) lead. 

FIG. 2: (color online) The spin polarization An^ at the left dot (a), Ariij at the right dot (b), 
nL + nR/2 (c), and the current J (d) as a function of the energy levels sl and cr in the two 
quantum dots. The parameters are: F^ = F^^ = 0.3, tc = T = 0.1, Uin = 20, Ugx = 5, and V = 1. 

FIG. 4: The current J vs. the spin bias V for two inter-dot interactions Uex = 5 (a) and Uex = 
(b). (c) The current J vs. the level €r and (d) the current J vs. the temperature T. The solid, 
dashed, and dotted curves are for levels at = and cr = —20, el = 0.5 and cr = —19.5, and 
cl = 1 and cr = —19, respectively. The other parameters are the same as in Fig.2. 



FIG. 5: The induced charge bias Ve (thin curves) and Vg/V (thick curves) vs. the spin bias V for 
intcr-dot interaction Uex = 5 (a) and Uex = (b). (c) Ve (i.e. Ve/V) vs. the level cr and (d) Ve 
(i.e. Ve/V) vs. the temperature T. The solid, dashed, and dotted curves are for levels at ex, = 
and cr = —20, cl = 0.5 and sr = —19.5, and e^, = 1 and cr = —19, respectively. The other 
parameters are the same as in Fig.2 



FIG. 3: (color online) (a) Schematic stability diagram of the DQD under the finite spin bias 
V. The thin dashed lines are the stability diagram while V = 0. (b), (c), and (d) are location 
schematics of energy levels in the (0, 1) , (t, 1), and (0, [) domains, respectively. 
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